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Graphene superlattices formed by periodic strain are considered theoretically. It is shown that 
electron energy spectrum consists of minibands obtained by folding of the cone at the boundaries of 
the superlattice Brillouin zone with very narrow anticrossing regions. Light absorption under direct 
interminiband transitions is shown to be frequency dependent and dichroic. Giant dichroism of 
absorption is demonstrated for doped graphene superlattices. Asymmetrical graphene superlattices 
act as quantum ratchets allowing for generation of photocurrent at absorption of normally incident 
light. The helicity-dependent photocurrent spectrum is calculated for doped superlattices with 
various asymmetries. 



I. INTRODUCTION 



II. ELECTRONIC STATES IN GRAPHENE 
SUPERLATTICES 



Graphene attracts a great deal of attention due to 
both interesting physical phenomena and device appli- 
cations. Physics of graphene possesses many intriguing 
effects in optics, transport and other fields^— An exam- 
ple is optical absorbance of graphene, which is a fun- 
damental constant in a very wide frequency range. Af- 
ter focusing on pure graphene, the interest is shifted to 
graphcne-based systems whose electronic properties can 
be changed by external parameters. One of the promis- 
ing avenues is strain engineering of electronic bandstruc- 
turei^ For example, gaps in energy spectrum open in 
strained graphene.^ In graphene, the strain is equiva- 
lent to gauge, or pseudomagnetic fields,— ^ which suppress 
anomalous magnetoresistance caused by weak localiza- 
tion^ 



Deformation of graphene results in periodic ripples 
with height ~ 0.1 ... 1 nm and a period ~ 10 nni i 10 ' 11 
Such periodically strained graphene samples have D2/1 
point symmetry which allows for dichroism of light 
absorption. Asymmetrical strain lowers symmetry of 
graphene to Ci v with the C2 axis along the strain di- 
rection. In such a system a ratchet effect is possible 
which means a generation of electric current under light 
absorption of unbiased sample. In particular, a contribu- 
tion to the photocurrent sensitive to the helicity of light 
normally incident on a strained graphene structure is al- 
lowed. Up to now such circular ratchet effect has been 
detected only in small graphene samples owing to scat- 
tering from sample edges ^ or for oblique incidence due 
to the presence of a substrate*^ 



In this work we investigate electronic states and calcu- 
late dichroic absorption in periodically strained graphene 
samples. For asymmetrical periodic strain, we study 
the circular ratchet effect and calculate the circular pho- 
tocurrent for various asymmetries of the strain. 



The Hamiltonian for strained graphene is given by^ 

H(p) = v a-(p±A), (1) 

where p is the quasimomentum operator, er = (a x , a y ) is 
a vector of Pauli matrices, and the two signs correspond 
to two valleys. We consider a graphene sheet periodically 
strained in the y direction: 

A x {y + d) = A x {y), 

where d is the period. The value of A x is modulated 
with the amplitude jh/d with the estimates^ yielding 
7 ~ 0.03 . . . 0.3. Owing to the periodicity, a superlattice 
(SL) is formed, and the envelope function satisfies the 
Bloch theorem 

*(x,y + d) = <f(x,y)c ikd , 

where k is the wavevector for motion along the y direc- 
tion. 

We consider the Kronig-Pcnncy model where strain is 
a piecewisc-constant periodic function of y. The wave- 
functions within the regions with the fixed value of the 
strain A x = Jih/d can be cast as 
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Here q is the wavevector for free motion along the x 
direction, fe, are determined by the energy E via fc, = 
^(E/hv ) 2 -(q- ~ii/d) 2 , and tan<^ 4 = ±ki/(q - Ji/d), 
where the upper and lower signs should be taken for con- 
duction and valence band states, respectively. The dis- 
persion equation has the form 



cos fed = (T u +T 22 )/2, 



(3) 



where T is the transfer matrix. For the periodic structure 
with TV layers in each period it equals to a product T = 
Tjv_iTjv_2 . . .TiTq, where Tj is the transfer matrix for 
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FIG. 1: Energy spectrum in a SL with a period formed by two layers with widths h = 0.3d, h = 0.7d, and the strain amplitudes 
71 = 0.7, 72 = —0.3. Three-dimensional plot of the spectrum (a), energy dispersion as a function of q for four selected values 
of k (b), dispersion along k for q = (c) and q = 1/d (d). 



ith strained graphene layer of the width U. It is denned 
as ^i{x, y + k) = T^^x, y) and given by 
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For a graphene SL with a period formed by two layers, 
the dispersion equation has the forrr^ 



COskd = COS fci^i COS &2^2 

+ sin kill sin &2^2 



(5) 



COS ifi COS <y?2 — 1 

sin</3i sin <y92 



Energy spectrum represents a result of folding of Dirac 
cone within one-dimensional Brillouin zone of the SL. 
The cone is shifted from q = by a value go = (71^1 + 
^2^2) I d 2 . The spectrum is given by conduction mini- 
bands E n (k,q) and valence minibands with the inverted 
dispersion —E n (k, q) with n = 0, 1, 2 . . . The conduction- 
and valence-band cones touch each other at q = qo, k = 0, 
i.e. the zero-energy state is degenerate. At the points 
k = 0, ±7t/g?, the gaps between the minibands are opened. 
At small deformation I71 — 72I <C 1, the gaps are given 
by 
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At arbitrary values of the strain in the layers, 7,, the 
energy spectrum is found from numerical solution of the 



dispersion equation ([3]). Figure [T] presents the energy 
spectrum in the SL with the layer widths li = 0.3d, I2 = 
0.7d, and the strain amplitudes 71 = 0.7, 72 = —0.3. One 
can see folding of the energy spectrum at the boundaries 
of the Brillouin zone of the SL (k = ±7r/rf), and the 
conical dispersion in the direction perpendicular to the 
strain. Figure [T] demonstrates that even the large strain 
with 7 ~ 1 results just in small deviations from the pure 
conical dispersion which occur in close proximity of the 
anticrossing points. Therefore the energy spectrum and 
the electronic states can be described in the nearly-free 
electron model where the strain is considered as a small 
perturbation. In this approach, the unperturbed energy 
spectrum in the nth conduction miniband is simply given 
by 



(7) 



where k n = k — (27rn/c?)sign{fc}. In the following for cal- 
culation of absorption in periodically strained graphene 
we will use the perturbation in the form 



= — «» V d 



assuming 7 <C 1. The energies and the wavefunctions in 
the nth miniband are given by 



E n = 25< 0) + £ 
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Here ^4°^ are unperturbed envelopes 

where the upper and lower signs are taken for the con- 
duction and valence minibands, and tan(^„ = k n /q. 
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III. DICHROIC ABSORPTION IN GRAPHENE 
SUPERLATTICES 

The Hamiltonian of electron-photon interaction in 
graphenc reads 



He-phot — -f • A, 

c 



(10) 



where A is the light-wave vector potential. This yields 
a general expression for absorbance under direct optical 
transitions from an initial (i) to the final (/) miniband: 

Vfi = 4 (2 ^° )2 \ e ■ "a I' S l E f (*,?)- Ei (fe, ?) - M • 

(11) 

Here the factor 4 accounts for the spin and valley de- 
generacies, u is the light frequency and e is the light 
polarization vector. 

Below in this Section we calculate the absorbance for 
interband transitions in undoped graphene SLs and for 
transitions between the ground and the first excited con- 
duction minibands for doped SLs. 



A. Interband absorption in undoped SLs 

We consider light absorption due to interband transi- 
tions between the ground valence and conduction mini- 
bands which have opposite dispersions: E v o(k,q) = 
—E c o(k,q). In the absence of strain, the interband 
absorbance has a universal value r\ = ire 2 /(hc). For 
normally-incident light, the Dih point group of the 
symmetrically-strained graphene sheet implies the follow- 
ing polarization-dependence of the absorbance: 



»7cO,«oH = T -{l-7 2 [$o(^) + *iHcos2a]}. (12) 



Here a is the angle between the light polarization vec- 
tor e and x axis. We take into account in Eq. © only 
the admixture to the cO and vQ states from the closest 
minibands cl and vl, respectively. Taking into account 
the strain both in the energies Eq. ([5J and in the optical 
matrix elements, we obtain the absorbance in the form 
of Eq. flT2J) with 
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FIG. 2: Polarization-independent ($o) and polarization- 
dependent ($1) absorbance spectra under interband transi- 
tions in the symmetrical undoped SL. Inset shows the per- 
turbation U (y) caused by periodical strain and the scheme of 
optical transitions. 
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Here r((p) = y/T+ 16C 2 — 8Csiny>, 
cos 2ip + 4£ sin <p 
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Spectral dependences of polarization-independent and 
polarization-dependent parts of the absorbance, $o(w) 
and $i(w), are presented in Fig. [2] One can see that 
the values of both parts have an order of 10 _2 7 2 . This 
means that, in samples with a large enough strain 7 ~ 1 5 
such modulation of absorbance can be detected experi- 
mentally. 

The both absorption parts, $o(w) and $i(w), should 
tend to constants at frequences uj — and uj = irvo/d 
due to anticrossing of the energy dispersions. We do 
not obtain this saturation due to oversimplified approach 
Eqs. ([S]). However the calculated spectra are correct 
everywhere beyond the narrow spectral regions Auj ~ 
jnvo/d -C 1 near these two points. 



B. Inter-miniband absorption in doped SLs 

Let us consider absorption caused by transitions be- 
tween the ground (cO) and the excited (cl) conduc- 
tion minibands. The polarization dependence of the ab- 
sorbance 77ci,c0 can be conveniently written as 



Va.co = Vx\e x \ 2 +Vy\eyf 



(16) 
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where e x , e y are the light polarization vector compo- 
nents perpendicular and along the strain direction, re- 
spectively. 

The absorbance is given by 

Th,v = 4^^2! £ |(clK,|cO)| 2 (/ c0 - f cl ) (17) 
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where f c o and f c ± are the occupations of the initial and 
final states. Here we put the unperturbed energies Eq. ([?]) 
because the matrix element of inter-miniband transitions 
already contains 7. The calculation yields 
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(2C 2 - l/2)/(l + 2Csin^), C is given by 



G x = sin 2 (<p + ip') , G y = 4 cos 4 

and the angle ip' is defined as follows 
, cosi^(s + l) 
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At low temperatures electrons are degenerate, and the 
light absorption under direct inter-miniband transitions 
is possible at E\ > Ep > Eq only, where Ep is the Fermi 
energy, see inset in Fig. [3] Therefore the occupation fac- 
tor is given by (/ c o— fd) — 0(EF/hu} — s)6(s+l — EF/huj), 
where 6{x) is the Heaviside function. 

The absorbance spectra r] y (uj) and rj x (uj) are presented 
in Fig. [3J One can see large dichroism of inter-miniband 



absorption: r\ x <C f)y. This is explained by different char- 
acter of electron motion parallel and perpendicular to 
the strain direction y. The light polarized along x can 
be weakly absorbed because of ballistic electron propa- 
gation in this direction. On the other hand, the inter- 
faces between strained and unstrained stripes play the 
role of coherent scatterers, allowing for absorption of the 
y-polarizcd radiation. High absorption dichroism demon- 
strated here is similar to the result of Ref. UH where grids 
of graphene nanoribbons were studied. In contrast, we 
show that even pure graphene samples can have highly 
dichroic absorption if they are periodically strained. 

The absorbance T] x is nonzero in graphene SLs in con- 
trast to conventional semiconductor SLs. The differ- 
ence is in the electron-photon interaction. According 
to Eq. ([TU|. the corresponding Hamiltonian in graphene 
contains the operator a x while in semiconductor SLs 
He-phot is proportional to the unit operator with zero 
matrix elements between the states in different mini- 
bands. 

Let us compare T] x ,y calculated above with the Drudc 
absorbance accompanied by scattering with the charac- 
teristic time t, rjo ~ (Ep /huj 2 T)(e 2 /he). For realistic 
parameters E-p ~ hu ~ 0.1 eV, t ~ 1 ps we have 
r)D ~ 10~ 2 e 2 /(fi.c), which is about an order of magni- 
tude smaller than rj y and comparable with r\ x at 7 ~ 1. 
This demonstrates that the absorption dichroism caused 
by the periodic strain is observable in real samples where 
the Drude absorption is also present. 



IV. CIRCULAR RATCHET PHOTOCURRENT 

Consider a graphene SL formed by asymmetrical strain 
lacking the reflection operation y — > —y, see insets in 
Fig. SJ This system has Ci v symmetry with the C2 axis 
along the y direction. Symmetry considerations imply 
that absorption of normally-incident circularly polarized 
light is accompanied by generation of a photocurrent 
flowing perpendicular to the strain direction. This means 
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FIG. 4: Circular-polarization sensitive photocurrent in asymmetrical SLs. Insets show strain profiles for different SLs. The 
difference between the highest and lowest strain in the layers, |7 ma:c — jmin] = 0.8, 0.45, 0.5, 0.4 and 0.3 for the structures 1. . . 5, 
respectively. The upper SL in the left panel coincides with the structure #3 in the right panel. The period d — 10 nm, the 
relaxation time in the cO miniband to = 1 ps, the Fermi energy _Ep = 3>hvo/d. 



that asymmetrical graphene SLs serve as quantum ratch- 
ets. 

The phenomcnological expression for the photocurrent 
density has the form 



perturbation of the form 



Jj 



XlPc: 



(19) 



where P c irc = i( e x e y — e*e y ) is the light circular polariza- 
tion degree. Microscopically, the density of the photocur- 
rent generated under inter-miniband optical transitions 
in doped graphene is given by 



ae^wffco (^n-^ro). (20) 



Here v enumerates the two valleys, c o ^ s ^ ne m ter- 

miniband optical transition rate in the ^th valley, Hv x n ^ = 
dE n /dq, and t„ is the momentum relaxation time in the 
nth miniband. The expression for x nas the following 
form: 

X = S -^J2{v^ 1 -vi% )^\M^M^} 
uj z hc kqv i J 

x(/ c0 - /ci)^i(fc,g) - E (k,q) - foj], (21) 

where M x "y = (ucl\a Xt y\vc0). 

Numerical calculations performed according to 
Eq. (f2~Tj) for various asymmetrical SLs show that under 
the change of the strain in all layers 7; — > 07^ the circular 
photocurrent scales as 



X -> a x- 



(22) 



This result can be confirmed by the analysis in the free- 
electron approximation with a non-centrosymmetric pe- 
riodic perturbation. For this purpose, let us consider the 



A x = 



71 cos ■ 



2ny ( 4iry \ 8iry 

+ 72 cos I —j- + p I + 74 COS —j- 



where /3 ^ 0,ir. It turns out that, for such perturbation, 
nonvanishing circular photocurrent is obtained only if all 
the three amplitudes, 71,72,74, are nonzero. Indeed, for 
direct optical transitions between cO and cl minibands, 
the transferred momentum 27r/d is taken up by the pe- 
riodic strain, so the matrix elements M x "y oc 71. In or- 
der to take into account the SL asymmetry, one should 
consider the 72-related contribution. However, the mo- 
mentum transfer is twice larger in this case, and the cor- 
rection to the matrix element can be 5M x "y oc 7172 only. 
The interference of these two corrections yield the 7^72 

contribution to Im | M x v ^ j which results in no elec- 

tric current after summation over two valleys since they 
differ by signs of 7's. Moreover, such 7-odd photocurrent 
in each valley would mean the valley current^ at circu- 
larly polarized excitation which is forbidden by symme- 
try. Therefore such a contribution nullifies already after 
summation over q and k. 

The only possibility to get a nonzero photocurrent is 
the fourth order in the strain. It can be obtained from in- 
terference of the first-order term M x v ) y oc 71 with a third- 
order contribution. However it is impossible to construct 
the correction to the cO — > cl transition matrix element 
which oc 7^72- The fourth-order in strain probability of 
the optical transition is obtained as a product of M x "y 
and 8M x u l oc 717274- This yields the circular ratchet 
current Eq. (|19p with x ~ 7i7274- This dependence co- 
incides with the numerical result Eq. (f22j). 

Figure U presents results of numerical calculation of 
the circular ratchet photocurrent for various asymmetri- 
cal SLs. We assume that the momentum relaxation time 
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in the 1st miniband is much shorter than in the ground 
miniband, Ti <C t , and disregard the corresponding con- 
tribution to the photocurrent. In calculation we take the 
period d = 10 nm, the relaxation time tq = 1 ps, and the 
Fermi energy Ep = 3hvo/d. 

Left panel of Fig. 0] represents the circular ratchet cur- 
rent for two strains with opposite asymmetry, i.e. for two 
SLs transforming one to another by a reflection y — > — y. 
The calculations demonstrate that the circular photocur- 
rents are directed oppositely and have the same magni- 
tude. The right panel of Fig. @] shows the circular pho- 
tocurrents for different strained SLs. One can see a larger 
photocurrent for SLs with higher degree of asymmetry 
For the symmetrical SL (structure # 5) the circular pho- 
tocurrent is zero. 



V. CONCLUSION 

High dichroic absorption in graphene SLs demon- 
strated above can be used in polarizers operating in in- 
frared and terahertz ranges. The magnitude of the cir- 
cular photocurrent is x ~ 1 n A cm/W at typical photon 
energy fouj ~ hvo/d ~ 60 meV for the period d = 10 nm. 



These estimates are close to experimental conditions re- 
alized in Refs. [Ull which allows one to expect obser- 
vation of the considered effects in experiments. 

To summarize, graphene SLs formed by periodic strain 
are considered theoretically. We show that electron en- 
ergy spectrum consists of minibands obtained by fold- 
ing of the cone at the boundaries of the SL Brillouin 
zone with very narrow anticrossing regions. Light ab- 
sorption under direct interminiband transitions is shown 
to be frequency dependent and dichroic. Giant dichro- 
ism of absorption is demonstrated for doped graphene 
SLs. Asymmetrical graphene SLs are shown to act as 
quantum ratchets allowing for generation of photocurrent 
at absorption of normally incident light. The helicity- 
dependent photocurrent spectrum is calculated for doped 
SLs with various asymmetries. 
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